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Abstract. The long-term variability of a sample of 149 
optically selected QSOs in the field of the Selected Area 
94 has been studied in the R-band. The relations be- 
tween variability and luminosity and between variability 
and redshift have been investigated by means of "robust" 
statistical estimators, allowing to disentangle the effects 
of the measurement errors. The results are compared with 
the corresponding properties of the variability in the B- 
band for the same sample. An anti-correlation between the 
R-band variability and the intrinsic luminosity is found, 
analogously to what is observed in the B-band. The am- 
plitude of the R-band variability turns out to be smaller 
(of a factor 1.13 ± 0.05) than the B-band variability. The 
implications in terms of the black-hole, starburst and mi- 
crolensing models are discussed. 

Key words: quasars: general - quasars 



1. Introduction 

The detailed study of statistically well-defined samples of 
optically selected QSOs is beginning to pin down their en- 
semble variability properties. An anti-correlation between 
variability and luminosity has been observed, such that 
more lumi nous QSOs show a smaller amplitude of the 
variations ( Hook 1994 , Cristiani et al. 1996). It starts also 
to emerge from the statistical noise a positive correlation 
between the amplitud e of the variability a nd the redshift 
( |Giallongo et al. 1991 L Trevese et al. 1994 , Cristiani et al. 
1996, Cid Fernandes et al. 1996 ). This latter dependence 
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may be interpreted as an anti-correlation between ampli- 
tude of the variability and w avelength in the individual 
objects (Cristiani et al. 1996, Cid Fernandes et al. 1996 , 
Di Clcmentc et al. 1996 ) , as suggested by observations of 



relatively small samp les of objects ( Perola et al. 1982 



Edelson et 
Courvoisier 
plications for 
generation 



1994) 



al. |199C| , Kinney et al. |l99l| , Paltani fc 



Both properties bear important im- 
thc understanding of the AGN energy 
mechanism and distinguishing among the 
competing variability models based on the black-hole 
(Rees 1984), starburst (Terlevich et al. 1992) and mi- 



crolensing (Hawkins 1993, 1996) scenarios. At the same 
time they are central for the quantification of the biases 
affecting samples selected on the sole basis of variabil- 
ity ( Veron & Hawkins 1995) or on the basis of multicolor 
data, if the photometry was carried out at different epochs 



in the various wavebands (e.g. Warren et al. 1991). 

However the variability-redshift correlation could be 
due to an evolutionary effect, in the sense of a higher ac- 
tivity in the early phases of the AGN life. The only way to 
remove such an uncertainty is to check directly the vari- 
ability properties of individual objects, by observing them 
simultaneously in two or more wavebands. For this reason 
we have investigated the R-band variability of the SA94 
QSO sam ple, the same that in Cristiani et al. 1996 (in the 
following |Paper 1 ) had been studied in the B band. 

We report in Section 2 the description of the SA94 
QSO sample and the photometric material used to inves- 
tigate its variability; in Section 3 the ensemble structure 
function is introduced as a statistical measure of the vari- 
ability properties and the detected trends are further anal- 
ysed with the aid of other statistical indices; the results 
and implications are discussed in Section 4. 

Throughout the paper we have assumed cosmological 
constants H Q = 50 Km sec -1 Mpc -1 , q = 0.5. 



2. The quasars sample 



2 
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2.1. Definition of the sample 

The SA94 sample is made up of 149 optically selected 
quasars S listed in Table 0, observed in a rectangular 
area of the sky within the limits 2 ,l 43 m 51.2 s < a < 
2 /l 59 m 14.3 s , -2°03'23.8" < S < 2°32'11.0" (Epoch 
1950.0), covering 17.66 sq. deg. 

They have been mainly discovered in the UVx and 



Table 1. List of the QSOs in the SA94 



objective-prism surveys described in Paper I and in the 
references to Table [|. 

2. 2. Calibration of the photographic material and error es- 
timation 

8 plates taken with the ESO La Silla and UK Schmidt tele- 
scopes and a copy of a POSS E-plate (E1453) have been 
analysed. In Table || a detailed list is given. The combi- 
nation IIIa-F + RG630 filter was used in all cases except 
for the plate E1453, that is a combination of a IIIa-E + 
a red Plexiglas 2444, providing a transmission similar to 
that of a Wratten n°29 filter (peak sensitivity near A6563). 
The natural photometric system defined by the IIIa-F + 
RG630 combination (in the following defined R') is close 
to the Johnson-Kron-Cousins R and the color transforma- 
tion is given by Bessel (1986): 



Ft! = R - 0.013 - 0.201 (R 
+0.0295 (R - I) 3 



I) + 0.100 [R-If + 



(1) 



The combination IIIa-E + red Plexiglas 2444 provides a 
system for our purposes indistinguishable from the pre- 
vious one, with no significant color term in the range 
0.5 < B — R < 3, as verified on our photometric mate- 
rial. 

The plate material has been scanned with the COS- 
MOS microdensitometer (MacGillivray & Stobic 1984) in 
IAM mode. The resulting tables, one per each plate, con- 
taining the instrumental magnitudes and other useful pa- 
rameters for the objects detected, have been merged to- 
gether in one table. Only objects with at least 3 detections 
in the 9 plates have been accepted in this final table con- 
taining 82 984 entries. The astrometric error box defining 
a common detection has been defined as a circle of 3.5 
arcsec of radius. In this way spurious detections (plate 
defects) are minimized to an acceptable level, while real 
measurements are in practice never discarded. 

The calibration of the magnitudes of each plate has 
bee n carried out according to the following procedure (see 
also Paper I ): 

1. the magnitudes of 108 standard stars, transformed into 
the natural colour system of the plates, R' , according 



1 In the present paper quasars are denned as objects with 
a starlike nucleus, broad emission lines, brighter than Mb = 
—23 mag, applying the K-corrections computed on the basis 
of the composite spectrum of Cristiani and Vio (1990), and 
dereddened for galactic extinction according to Burstein and 
Heiles (1982). 
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2.521 


17.48 


+0.00 


a 


2 


58 


11.37 


+0 


05 


06.6 


1.727 


18.63 


-0.01 


d 


2 


58 


11.53 


+0 


09 


42.8 


1.497 


19.35 


+0.23 


d 


2 


58 


14.54 


+0 


42 


50.7 


0.661 


18.46 


+0.06 


d 


2 


58 


14.72 


+1 


37 


06.2 


1.302 


18.74 


+0.04 


d 


2 


58 


25.73 


+1 


37 


39.1 


0.595 


18.63 


+0.00 


d 


2 


58 


54.46 


+1 


45 


50.4 


1.349 


18.87 


+0.04 


d 


2 


59 


02.77 


+1 


12 


53.7 


2.316 


18.74 


+0.00 


d 


2 


59 


03.14 


+1 


26 


27.3 


1.578 


18.38 


+0.06 


d 


2 


59 


06.29 


+1 


04 


03.8 


1.770 


18.48 


+0.01 


d 
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Table 2. List of the Plates 



Plate 


Exp 


Date 


Epoch/ date 


Limit 




(min) 






(mag) 


E1453 


40-60 


1955 Oct 22 


1 (1955.81) 


19.5 


R5413 


60 


1983 Dec 08 


2 (1983.94) 


20.3 


R5414 


60 


1983 Dec 08 


2 (1983.94) 


20.3 


R12218 


150 


1987 Oct 16 


3 (1987.79) 


20.1 


R12807 


140 


1988 Oct 15 


4 (1988.79) 


19.5 


R13409 


100 


1989 Oct 28 


5 (1989.83) 


20.1 


R13420 


100 


1989 Nov 01 


5 (1989.83) 


20.1 


R8757 


120 


1990 Sep 21 


6 (1990.72) 


19.1 


R9477 


120 


1991 Oct 07 


7 (1991.77) 


17.3 
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Fig. 1. QSO rest-frame structure function in the R-band Fig. 2. Structure function for QSOs in the R-band (continuous 

line) and in the B-band (dashed line) 



to Eq. 1, have been used to derive for each plate a poly- 
nomial regression between instrumental and calibrated 
magnitudes. 

2. the median magnitude of all the objects has been com- 
puted and each plate has been re-calibrated against the 
median magnitudes of the 27111 point-like objects. 

3. a procedure of uniformization of the usually spatially 
variable response of the photographic plates has been 
applied. Each plate has been subdivided in 10 x 10 sub- 
areas; for each of them the differences between the ref- 
erence and the individual plate magnitudes have been 
computed and their distribution analysed. The zero- 
point shifts estimated in this way for each sub-area as 
a function of the magnitude have been smoothed and 
applied to the re-calibrated magnitudes. In the follow- 
ing we will refer to the magnitudes obtained in this 
wayasify w . 

4. for each plate the uncertainties of the R'f ina i magni- 
tudes have been estimated by analysing the distribu- 



tion of the differences AR = R 



final 



R', 



median 



as a 



function of the median magnitudes for all the point- 
like objects. 

5. the completeness limit of each plate (see Table 2) has 
been estimated as the maximum of the histogram of 
the calibrated magnitudes. 



3. Statistical indices of the variability 



3.1. The structure function 

To study the ensemble variability properties of the 
SA94 sample we have used the structure function 
( [gimonctti ct al. 1985 , Paper I), defined as: 

SF{t) = {[mag(t) - mag(t + r)] 2 ) (2) 

where mag(t) are the magnitudes at the time t, r is a 
time lag evaluated in the rest-frame of each QSO and the 
brackets " (...) " indicate a mean over the ensemble. The 
effect of the measurement errors has been subtracted ac- 



cording to the procedure described in detail in Paper I 



Only objects detected in at least 5 epochs and with a me- 
dian magnitude brighter than the completeness limit have 
been used to compute the SF. Only 6 QSOs have been 
measured on Plate R9477 brighter than the completeness 
limit. We have checked that the SF and the following sta- 
tistical results do not depend on the inclusion/exclusion 
of these data. 

In Fig. [I] the R-band SF in the quasar rest frame is 
shown. We observe a steady rise, steeper up to about 
5 years and a flattening afterwards (within the errors), 
reaching values corresponding to about 0.1 mag 2 of RMS 
variability. 

It should be noted that points in Fig. |l| with At > lOyr 
are all derived from the comparison of the POSS plate 
with the remaining plates. Being 28 years apart from Plate 
R5413, Plate E1453 does not play any role in the estimates 
of the SF on the shorter timescales (5 yr or less). 

In Fig. H a comparison between the SFs in the B-band 
(as derived from the data discussed in Paper I) and in the 
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Fig. 3. ABi vs. ARi (see text). Dashed line: straight line with 
a slope 1. Continuous line: least-squares fit between ABi and 
ARi (see text). The cross in the lower-right corner shows the 
typical measurement errors 



R-band is reported. The SF in the B-band shows a larger 
variability amplitude with respect to the one in the R- 
band. To quantify the statistical difference between the 
two SFs, we have calculated the weighted average of the 
ratio SF K /SF B , for lags between 1 and 5 yr in the QSOs 
rest frame. The result is SF R /SF B = 0.66 ±0.13, a differ- 
ence from the unity significant at a 2.5a level. 

In order avoid any possible bias deriving from differ- 
ent time samplings or time dilation effects, a direct com- 
parison of the variability in the R and B band has been 
carried out. As shown in Fig. |^, for each object we have 
computed the quantities ARi — Ri — (R) , as the difference 
between the magnitude at a given epoch and the average 
magnitude. Since the epochs of the R plates are in large 
part coincident with the epochs sampled with the B plates 
(epochs 2, 3, 4, 5, 6, 7 of the present paper correspond to 
epochs 2, 6, 8, 9, 10, 11+12 of Paper I, respectively) it 
has been possible to plot the ABi versus the ARi for all 
the objects with at least 5 epoch-magnitudes in common. 
Again, the POSS plate does not play any role in this com- 
parison. 

We have then fitted a straight-line relationship be- 
tween ABi and ARi, with a least-squares technique that 
takes into account the errors of the data points on both 
axes ( |Fasano fc Vio 1988[ ). The result is AB = (1.130 ± 



Fig. 4. Distribution of the SA94 QSO sample in the M B — z 
plane 



ditional assumption that the R flux varies in phase with 
the B flux, i.e. that for AR — AB — 0, then we obtain 
AB = (1.130 ± 0.054)Ai?, again fully consistent with no 
intrinsic scatter. The coefficient between the R and the B 
variability obtained in this way is consistent within the 
errors with what has been estimated on the basis of the 
SF analysis. 

In Fig. |i| the distribution of QSOs in the redshift- 
magnitudc plane is illustrated. To further investigate and 
disentangle the dependences of the variability on the lu- 
minosity and redshift we have examined two subsamplcs. 
The first is defined within the redshift limits 0.3 < z < 1.4 
(the area between the dashed lines in Fig. ||). The second 
is defined within the absolute magnitude limits —26.85 < 
Mb < —24.50 (the area between the continuous lines in 
Fig.0). 



As in Paper 1 



0.054)Ai? + (0.004 ± 0.009). The dispersion of the points 
around this relation is fully compatible with the measure- 
ment errors, with no intrinsic scatter. If we make the ad- 



the SFs indicate that the QSO 
variability (in magnitude) is, also in the R band, clearly 
anticorrelated with absolute luminosity, in the sense that 
more luminous objects have a smaller variability ampli- 
tude (Fig. |). No correlation between redshift and vari- 
ability - of the type discovered in Paper I for the B-band 
variability - has been detected within the errors (Fig. ||). 
However, such a failure is not surprising considering that 
in the present work we are dealing with only one sample 
(in place of the three of Paper I) and the still non-uniform 
coverage of the redshift-absolute magnitude plane of the 
subsample defined between —26.85 < Mb < —24.50 in- 
troduces a luminosity-redshift correlation in the data that 
may mask any redshift-variability correlation. 
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Fig. 5. The structure function for quasars of different lumi- 
nosity in the redshift interval 0.3 < z < 1.4. Continuous line: 
QSOs with M B < -24.5. Dashed line: QSOs with -24.5 < M B 



Fig. 6. The structure function for quasars of different redshifts 
in the luminosity interval —26.6 < Mb < —24.5. Dashed line: 
QSOs with z < 1.5. Continuous line: QSOs with 1.5 < z 



3.2. The variability index 

In order to further verify the results obtained in the previ- 
ous section, we have evaluated the variability index (IDX) 
for each QSO, defined as the variance of the normal pro- 
cess representing the intrinsic variability that, added to 
the photometric errors, reproduces the observed light vari- 



ations on a given range of timescales (see Paper I). The 
index has been evaluated in one time bin (from 1 to 15 
years in the QSOs rest frame) to avoid on the one hand 
time scales for which the structure function is still quickly 
rising, as observed in the previous section, and on the 
other hand time lags unaccessible with the present data 
for the higher-redshift QSOs. Then we have analysed the 
correlation matrix, reported in Table ul 



significance greater than 99%). Due to the strong flux- 
rcdshift anticorrelation in the QSO sample it is not pos- 
sible to disentangle if one of the two correlations is spuri- 
ous. However, if, on the basis of the results of the previous 
subsections and of Paper ] , we assume as fundamental the 
absolute magnitude - variability correlation, it is possible 
to subtract its influence from the variability-redshift anti- 
correlation via the method of partial correlation analysis 
( [Spiegel 1991 ) . Applying this recipe to the results of Ta- 
ble the corrected correlation coefficient between vari- 
ability index and redshift, turns out to be 0.03 ± 0.08 
(65%), not significant, in agreement with the result shown 
in Fig. |6[ 

4. Discussion 



Table 3. Correlation coefficients 





R 


z 


Mb 


IDX 


0.09 ±0.08 


-0.18 ±0.08 


0.24 ±0.08 


M B 


0.42 ±0.07 


-0.81 ±0.03 




z 


0.15 ±0.08 







The correlation coefficients indicate a correlation of 
the variability index with the absolute magnitude and an 
almost equal anticorrelation with redshift (correlation co- 
efficients respectively of ±0.24 and —0.18, both with a 



The decrease of the variability amplitude with increasing 
wavelength observed in the present QSO sample is con- 
sistent with the results obtained for individual objects at 



low redshift by Edelson et al . ( |1990| ), Kinney et al. ( |199lD , 
Paltani & Courvoisier (1994). Is this variability-frequency 
correlation sufficient to account for the po sitive cor relation 
between redshift and variability found in Paper I ? 

Let us assume a simple power-law representa- 



tion of a typical QSO spectrum ( Giallongo et al. 1991 
Pi Clcmcntc ct al. 19961) : 



Fx cx 



(3) 
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where Fx is the flux at an observed wavelength A, Ao an 
arbitrary wavelength and a the spectral index. The vari- 
ability can then be modeled as a change in the spectral 
index: 



a+e(t) 



(4) 



where e(t) is a function of time. In this way the flux varia- 
tions are idealized as changes in the slope of the power-law 
"hinged" on the point of wavelength Ao, that we imagine 
somewhere redwards of the B and R bands. 
Converting in magnitudes, 



Mr (*) = 2.5 log 



Ar 



A (1 + z) 



a+e(t) 



+ c, 



(5) 



we obtain the following expression for the magnitude vari- 
ations: 



AM R = M R (t 2 ) - M R (h) = 
= 2.5[e(i 1 )-e(t 2 )]log 



A (l + z) 



Ar 



(G) 



Applying the same equation to the B-band, we can write: 



A Mb _ log^ + log(l + ^) 
AM R ~ log ^ + log (1 + z) 



(7) 



From Eq. we would expect a decreasing ratio 
AMb/AMr with increasing redshift. In the present data 
we observe, if any, a very small decrement of AMb / AMr 
with the redshift (from 1.14 ± 0.08 to 1.10 ± 0.09 for two 
sub-samples with redshift z < 1 and 1 < z < 1.8 respec- 
tively) . 

Introducing the value of AAf B /AM R = 1.13 ± 0.05 at 
a (z) ~ 1.5 reported in section 3 in Eq. [?], allows us to put 
a lower limit to the wavelength Ao ^ 31 000 A. 

If we turn back to Eq. 6, we can see how this lower 
limit to Ao, together with the typical observed RMS for 
the variability AMb arid AMr of 0.35 and 0.28 mag, re- 
spectively (see Fig. ||), puts an upper limit to the RMS 
variation of the spectral index a e ^$ 0.11. 

This value can be compared with the typical disper- 
sion in the spectral indices observed in various QSO sam- 
ples: Pei et al. (1991) have obtained a a ~ 0.24 in a control 
sample of 15 QSOs, Francis (1993) observed a a ~ 0.5 in 
a much larger sample, showing also a systematic harden- 
ing of the spectrum with increasing redshift. In this way 
the dispersion in the observed spectral indices of the QSO 
population can be only partially accounted for by vari- 
ability and significant intrinsic differences between QSO 
and QSO have to be invoked (but see also the dispersion 
G a om ~ 0.1 estimated by Elvis et al., 1994). 

Putting in Eq. 6 our best estimates for er e = 0.07 and 
Aq = 86 000 A, we obtain for the typical RMS variability 



in the B band 0.26 and 0.33 mag at z = 0.5 and z = 3, 
respectively and 0.22 and 0.30 in the R band at the same 
redshifts, consistent with the observations of Paper I. 

The dependence of the variability time scale on the 
wavelength is a standard prediction of the accretion disk 
model. The temperature of the disk decreases with the ra- 
dius, and the dynamical, thermal and viscous time scales 



all increase as r 3 / 2 (see Baganoff & Malkan 1995 and ref- 
erences therein). In particular, in the standard model, 
involving instabilities in a thin accretion disk around a 
supermassive black hole, the intrinsic characteristic time 
scale should vary as r oc A 2 . More in general, if the source 
is thermal and if the flux variability is due to temper- 
ature changes, the observed wavelength depe ndenc e can 
be naturally produced (Paltani & Courvoisier 1994 ). The 
only requirement concerns the turnover in the spectral 
energy distribution that should not be located at such 
a high frequency that the spectral slope in the observed 
band becomes independent of the temperature (e.g. the 
Rayleigh- Jeans domain for a black body spectrum). As 
remarked by Di Clemente et al. (1996), if brighter objects 
are on average hotter and the variability is due to tem- 
perature changes, one would expect a negative correla- 
tion between amplitude of the variability and luminosity, 
as observed. Assuming for example a spectrum of black 
body or thermal bremsstrahlung, the spectral turnover of 
brighter objects is progressively shifted at higher frequen- 
cies producing progressively smaller flux changes, from 
SI/I oc 5T/T{hv/kT) for kT < hu, to SI/I cx ST /T for 
kT > hv. 

The broad-band variability in the starburst model 



( Axetxaga et al. 1996 and references therein) is defined by 
the superposition of a variable component, supernova ex- 
plosions (SNe) generating rapidly evolving compact super- 
nova remnants (cSNR), and a non- variable component, a 
young stellar cluster and the other stars of the galaxy. The 
problem of predicting the wavelength dependence of the 
variability can be reduced to computing the spectra of the 
variable and non-variable components, and their relative 
luminosities. The spectrum of the non- variable component 
has been predicted to show a F v oc v a with a ~ — 1.0 -f- 
( Cid Fernandes fc Terlevich 1995 ). The variable/non- 
variable relative luminosities can also be estimated on 



the basis of stellar evolution ( Aretxaga fc Terlevich 1994 ). 
Much more difficult is to predict the optical/UV spectrum 
of the variable component. It is expected to be harder than 
the non-variable component with a F v oc v a and a ~ 0.5. 
We can model the total QSO flux according to a Simple 
Poissonian model ( Cid Fernandes et al. 199q ) as 



F, 



F mT + F\ 



bck 



(1 - /bck) 



Ao 



(a w -2) 



+ /l 



bck 



Ar 



(«bck-2) 



(8) 



where /b c k is the non-variable fraction of the total flux 
and abek, aw are the slopes of the non- variable and vari- 
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able components. Indicating with Ab and Ar the effective 
wavelengths of the B and R-band and putting for simplic- 
ity Ao = Ab, we obtain for the ratio between the B-band 
and R-band variability: 



A Me 

Mr 



/bck (^) + /var(l + Z) 
/bck + /var(l + z) A " 



(9) 



with Aa = a bc k - a var- 

For example, with Aa = 1 and /b c k, /var = 0.7, 0.3 re- 
spectively ( Cid Fernandes et al. 1996 ), the ratio between 
the B-band and R-band variability is expected to be about 
1.4 at z = 0, tending to 1 for increasing redshift, indepen- 
dent on the intrinsic luminosity of the QSO. The above 
reported slight tendency for a decreasing AMb / AMr ra- 
tio with increasing redshift is compatible with this model. 
However, an increase of the AMb /AMr with the ab- 
solute luminosity is also observed. We find the values 
AMb/AMr = 1.08 ± 0.06 and AMb / AMr = 1.31 ± 0.11 
for two subsamples with Mb > —25.2 ((z) = 1.05) and 
M B < -25.2 ((z) = 1.92), respectively. This suggests that 
more general situations have to be envisaged in which the 
fraction on the non- variable component and/or the pulse 
properties (e.g. the spectral distribution) depend on the 
global luminosity of the object. 

Recently Hawkins ( 1996] ) has proposed that the vari- 
ability of nearly all the QSOs (except the very low redshift 
ones) is produced by microlensing. A non-achromatic be- 
haviour can be accommodated in this kind of model. For a 
source comparable in size to the Einstein ring of the lens, 
with a quasar disk redder at larger radii, the bluer com- 
pact core would produce a larger amplitude of the varia- 
tions in the blue passband whereas the larger extent of the 
red image would cause a small amplitude variation, albeit 
with a longer duration, in the red. However it cannot be 
denied that achromatic light variations would have been 
one of the most significant characteristics of the microlens- 
ing model and the present observation of the wavelength 
dependence make it less appealing/compelling. In partic- 
ular the superposition of the two competing effects of time 
dilation and correlation of the variability with frequency, 
i.e. redshift, accounts naturally, also in the accretion disk 
and in the starburst scenarios, for the observed lack of 
correlation in the observer's rest frame between the vari- 
ability time scale and the redshift, otherwise claimed as 
an evidence in favor of the microlensing model. 

Future detailed comparisons between simulated and 
observed light curves will be fundamental to disentangle 
among these three scenarios and to start exploring the 
parameter space of each model. 
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